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ABSTRACT 



Aims. The aim of this paper is to analyse the validity of the Dopplergram and A-meter techniques for the Doppler 
diagnostics of solar photospheric velocities using the Ball 4554 A line. 

Metliods. Both techniques are evaluated by means of NLTE radiative transfer calculations of the Ban 4554 A line in 
a three-dimensional hydrodynamical model of solar convection. We consider the cases of spatially unsmeared profiles 
and the profiles smeared to the resolution of ground-based observations. 

Results. We find that: (i) Speckle-reconstructed Dopplergram velocities reproduce the "true" velocities well at heights 
around 300 km, except for intergranular lanes with strong downflows where the velocity can be overestimated, (ii) 
The A-meter velocities give a good representation of the "true" velocities through the whole photosphere, both under 
the original and reduced spatial resolutions. The velocities derived from the inner wing of smeared Ban 4554 A line 
profiles are more reliable than those for the outer wing. Only under high spatial resolution does the inner wing velocities 
calculated in intergranular regions give an underestimate (or even a sign reversal) compared with the model velocities, 
(iii) NLTE effects should be taken into account in modelling the Ban 4554 A line profiles. Such efi'ects are more 
pronounced in intergranular regions. 

Conclusions. Our analysis supports the opinion that the Dopplergram technique applied to the Ban 4554 A line is a 
valuable tool for the Doppler diagnostics of the middle photosphere around 300 km. The A-meter technique applied 
to this line gives us a good opportunity to "trace" the non-thermal motions along the whole photosphere up to the 
temperature minimum and lower chromosphere. 
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First studies of the barium lines in stellar spectra started 
in the 19 30s. One of the examples o f these earlier stud- 
ies is by iBidelman fc Keenaiil (|l95ll ). who pointed out 
that the anomalous strength of the Ball 4554 A line in 
some stars could be explained by a deviation from lo- 
cal thermodynamic equilibrium (NLTE). At the present 
time stellar studies are aimed at the abundance de- 
termination of barium isotopes using one-dimensional 
(ID) plane- parallel atmospheric modej-S and the NLTE 
assumption (iGigasl 1988; Mash onkina fc Gehrenl. 2000[ 
Mashonkina et al.l. l2003„ .2008 : .Mashonkina fc ZhaoL l200e 



Short fc HauschildB . l2006f ). Such studies play important 
role in estimating the yields of s- and r-processes in the 
nucleosynthesis of heavy elements in the Galaxy. 

Barium lines in the solar spectrum have been stud- 
ied since 1960s. Until the mid-seventies most the stud- 
ies aimed at determining the solar barium abundance us- 
ing a ID+LTE appr oach (Goldberg, Miiller, fc Alleii [1960|; 
iLambert fc Warneil . [1968 : .Holweger fc MiiUen. ,19741). The 
only e xc eption was the publication o f iTandberg-HanssenI 
(fT96l : iTandberg-Hansse n fc Sm^^ (I1970D . who had 
shown the importance of NLTE effects for the forma- 
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tion of the Ban 4554 A line. Later, iRutte^ (|1977l . |1978[ ): 
iRutten fc MilkevI (|l979t ) tackled in detail the NLTE Ban 
line formation probl em. Empirical analyses of centre-to- 
limb observations by iRuttenI ([l97l) have shown that the 
source function of the Ban 4554 A line deviates signifi- 
cantly from the Planck function. The effects of partially 
coherent scattering have also to be taken into account in 
order to reproduce observations away from the disc cen- 
tre. The hyperfine structure and isotopic shift play a very 
important role in the analysis of this line as well. 

For several reasons, a new debate has arisen in the 
recent literature on the solar Ban 4554 line. Firstly, a 
more realistic representation of the solar atmosphere by 
three-dimensional (3D) hydrodynamical simulations has 
become available. Relaxing the constraints of the plane- 
pa rallel (ID) modelling a new LTE abundance analysis 
by lAsplund et"an (|2005h resulted in a photospheric bar- 
ium abundance of = 2.17, close to the meteoritic 
value = 2 - 16. W hile the recent NLTE analysis by 
lOlshevskv et al.l (|2008l ) based on a 3D model lowers the 
value of ^Ba strictly to the meteoritic one. Note that the 
classical ID-approach resulted in a rather wid e spread of 
the solar barium abundance = 2.09 (R oss fc Alleii 
[l976l). A n. = 2.13 dG revessel Il984l ). and ^Ba = 2.40 
(|Ruttenl . [T978h . 
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Secondl y, according to th e atlas of the "Second Solar 
spectrum" (jGandorferl . l2002f ) . the linear polarization Q/I 
of this line close to the limb (/x = 0.1) is very strong (0.6%). 
The recently published theoretical investigation on the role 
of resonance scattering and magnetic fields in the polariza- 
tion sig nals of both B a n 493 4 A, and 4554 A resonance 
lines by iBelluzzi et al.l (|2007f) has demonstrated their im- 
portance for the measurements of weak magnetic fields on 
Sun and stars. 

Finally, several properties of the Ban 4554 A line have 
drawn attention to it as a diagnostic tool for the velocity 
field of the solar atmosphere. Due to the large atomic weight 
of barium (137.4 a.u.), one might expect low sensitivity of 
the line opacity to temperature variations and line-width 
insensitivity to thermal broadening. In addition, the Ban 
4554 A line has steep wings and a deep core. In a stan- 
dard ID model the core of this line is formed around 700 
km in the chrom osphere whil e the wings are photospheric 
(lOlshevskv et al.L 2008; Siitt crhn et al.l . 120011 ). As a result, 
the Ban 4554 A line gives an excellent opportunity to 
"trace" non-thermal motions (granulation and supergranu- 
lation velocity field and waves) th roughout phot osphere and 
even in the lower chromosphere. iNoved ()1967[ ) was one of 
the first to draw attenti on to th e Dopp ler diagnostic poten- 
tial of this line. Later. iRuttenI ()1978 ) confirmed the Ban 
4554 A line to be a perfect tool for investigating the ve- 
locity structure of the solar photosphere and lower chromo- 
sphere. Recently. [Siitterlin et al.l ()200lD have made the first 
serious attempt to use the Ban 4554 A line for mapping 
the line-of-sight (LOS) velocities (Dopplershift map) of dif- 
ferent structures in the solar photosphere. They presented 
observations with the Dutch Open Telescope (DOT) testing 
the D opplergram capabili t y of n arrow-band (80 mA) Lyot 
filter (ISkomorov skv et"aLl . Il97 6^ imaging the solar surface 
in the wings of this line in combination with speckle recon- 
struction. The Ba n 4554 A line is found to be an excellent 
tool for high-resolution Doppler mapping. 

Summarizing all the above, the Ban 4554 A line pro- 
vides a valuable diagnostic tool for the solar and stellar 
atmospheres. It is thus extremely important to investigate 
carefully the validity of the different data-processing tech- 
niques and interpretations applied to this line to obtain in- 
formation on physical conditions in the solar atmosphere. 
Our paper presents an example of such an investigation. 
Below we analyse two techniques used for recovering the 
solar velocity field from observations in this line. The aim 
of our analysis is to consider the advantages and disad- 
vantages of these techniques and to evaluate to what ex- 
tent the LOS velocities provided by them give a correct 
measure of the solar values. We keep in mind that, be- 
sides the techniques themselves, there is another important 
source of uncertainties that can affect Doppler diagnostics. 
The ground-based observations are typically affected by the 
Earth's atmospheric turbulence (seeing) and instrumental 
effects due to light diffraction on the telescope aperture (the 
finite spatial resolution of the telescope), the finite instru- 
mental width of the filters used, stray light, etc. So we tackle 
the problem taking into account seeing and instrumental ef- 
fects. The first technique analy sed in this pape r is a 5 -point 
Dopplergram method used by ISiitterlin et all (|2001h . The 
second technique i s known as a A-meter, first proposed by 
IStebbins fc Good3 (|l987t) . In subseq uent years it has been 
adopted by several researchers (see iKostik &: Khomenkol 
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and more references therein) for spectral observations of 
different Fraunhofer lines. 

The organization of this paper is as follows. 
Subsections 2.1 and 2.2 describe the 3D hydrodynamical 
model, atomic data and numerical methods needed for 
the NLTE radiative transfer calculations with the barium 
atomic model. In Subsections 2.3 and 2.4 we discuss the 
procedure used to simulate spatial smearing of the two- 
dimensional maps of the synthetic Ban 4554 A line profiles. 
We focus on seeing and instrumental effects. We establish 
theoretical calibration dependences between the granula- 
tion contrast and the Fried parameter (which specifies the 
characteristic size of atmospheric turbulence cells) based 
on a 3D-approach. Subection 2.5 defines artificial datasets 
employed in the paper while Section 3 presents obseva- 
tions. In Section 4 we discuss the validity of the 5-point 
Dopplergram technique used to obtain LOS velocities from 
the speckle-reconstructed observations of the Ban 4554 A 
line. Section 5 presents the results for the A-meter method. 
We discuss the granulation velocity field that one would 
expect from observations of this line under perfect spatial 
resolution (Subsection 5.1) and under different seeing con- 
ditions of the ground-based observations (Subsection 5.2). 
An extra point of particular interest has been to estab- 
lish heights from which the information on the velocity and 
intensity variations originate. We discuss this problem in 
Subsection 5.3. Finally, Section 6 presents our conclusions, 
while the Appendix gives a brief description of our NLTE 
modelling with emphasis on the NLTE mechanisms of for- 
mation of the Ban 4554 A line. We show population depar- 
ture coefficients, NLTE source functions and profiles of this 
line to illustrate the difference in NLTE results for granules 
and intergranules. 

2. Method 

2.1. 3D atmospheric model 

We use a 3D snapshot from realistic radiation hydrody- 
nami c al simulations o f solar convect ion (jStein fc Nordluiidl . 
119981: lAsplund et all . Il999l l2000airbl ). This simulation is 
based on a realistic equation of state, opacities and de- 
tailed radiative transfer. The size of the simulation box is 
6 X 6 X 3.8 Mm, with 1.1 Mm being located above the con- 
tinuum optical depth equal to one. To reduce the amount of 
time-consuming NLTE radiative transfer calculations, the 
snapshot was interpolated from the original resolution of 
200 grid points in the horizontal direction to a coarser res- 
olution of 50 grid points. At the same time, we increased 
the resolution in the vertical direction, taking only the up- 
per 1.1 Mm part of the snapshot between —200 km and 
900 km and interpolating from 82 to 121 grid points. Thus, 
the final 3D model has 50 x 50 x 121 grid points or 50 x 50 
ID models. 

It was concluded from previous studies that the 
3D model used here performs very satisfactory in 
terms of spectral line formation both for lin e shape 



2004; 



2001 



and asymmetrie s (lAsplund et al.l. Il999l l2000allbl. 

Shchuki na et al.l. 20051: Shchukina fc Truiillo Buenq . 

I2009t iTruiillo Bueno et all . l2004f). The model reproduces 
all the ma in features of solar c onvec tive velocities and 
intensities (|Kostyk fc Shchukinal . |2004( ). It reproduces as 
well the root-mean-square (r.m.s.) contrast of solar gran- 
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Fig. 1. Grotrian diagram of an atomic model for for Bai 
+Baii. Solid lines indicate radiative bound-bound transi- 
tions. The Ban A4554 A line arises from the 6s^Si/2 — 
6p^P3/2 transition. 



ulation at the 6301 A continuum wavelength (as ob- 
tained by HINODE), centre-to- limb variation of contin- 
uum inte nsity, and the polarization of the solar continuum 
(|Truiillo Bueno fc Shchukinal . |2009| ). 



2.2. Spectral synthesis 



According to iRuttenI (|1978[ ): iRutten fc MilkevI (|1979D ; 
lOlshevskv et al.l (|2008[) NLTE is the crucial factor that 
should be taken into account when calculating the 
Ball 4554 A solar line profile. Here we considered 
the NLTE barium line formation problem neglecting 
the effects of horizontal radiative transfer (1.5D ap- 
proximation). A self-consistent solution of the kinetic 
and radiative transfer equations has been obtained with 
an effici ent multilevel transfer code "N AT A J A" devel- 
oped by IShchukina fc Truiillo Buenol (|2001| ) to facilitate 
NLTE radiative transfer simulations with very complex 
atomic models. Before the code was successfully used 
for NLTE interpretati on of iron, oxygen, titaniurn, and 
stron tium solar s pectra (IShchukina fc Truiillo Bu enol.l200l[. 

ini, 



20091: [Shchukina et all. 120 05; Truiillo Bueno fc Shchukin, 
2007t [T ruiillo Buen o et al.i . .2004: Kostvk et al... ,2006t 



Kostvk fc Shchukina 
code is based 
calculations 



20011) . The 



2004 [Khomen ko et all 
on iterative methods for radiative transfer 

1995i : 



(seelTruiillo Bueno fc Fabiani Be ndicho. 

ISocas-Navarro fc Truiillo Buenol . 19971 and more references 
therein) that allow a fast and accurate solution of NLTE 
transfer problems. 

Our atomic model includes 40 energy levels of Ba i and 
Ban (see Fig.[T]). Note that all levels with hyperfine struc- 
ture (HFS) were treated as a single level, i.e. for the HFS 
sublevels the population departure coefficients (3 were taken 
to be equal. The levels are interconnected by 99 bound- 
bound and 39 bound-free radiative transitions. All levels 
are coupled via collisions with electrons. The atomic model 
and atomic data, including oscillator strengths, bound-free 
cross-sections, electron collisional rates, etc., are described 
in detail byfc)lshcvsky ct al. (2008). 

The departure coefficients f3 found from the self- 
consistent solution of the kinetic and radiative transfer 
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Fig. 2. Hyperfine and isotopic structure of the Ba n 4554 A 
line. Top: splitting profile patterns (thin solid lines) and re- 
sulting profile (thick solid line) in the ID MACKKL model 
(jMaltbv et al.l . il986l ). Two blends are shown by dotted (Cri 
4553.945 A) and dashed (Zrn 4553.970 A) fines. Middle: 
equivalent widths and isotopic shift of the patterns. Bottom: 
bisector due to isotopic shift (small open circles) calculated 
in the MACKKL model and bisector of the spatially av- 
eraged profile (3D model) caused by the granular velocity 
field (sofid line). 



equations were used as input to carry out the formal so- 
lution of the radiative transfer equation for the Ba n A4554 
A line. At this step of the solution we took into account 
the hyperfine structure (HFS) and isotopic shift of the Ban 
A4554 A line (see Fig. [2]). Note that because of non-zero nu- 
clear spin only two isotopes (^'^^Ba and ^■^^Ba) have HFS- 
splitting. The energies of HFS-sublevels are calculated ac- 
cording to a formula given by Radzig and Smirnov (1985). 
The hyperfine structur e constant A n eeded for these calcu- 
lations was taken from lRuttenI (|l978l ). We neglected the in- 
teraction of electrons with a nuclear electric quadrupole mo- 
mentum of because this effect is rather small. The isotopic 
shift was derived using mass shift and field shift constants 
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from lBerengut et al.l (|200"3l) and m ean square nuclear radii 
from lSakakihara fc Tanakal(|2001h . We synthesized the Ban 
A45 54 A line profiles employin g the isotopic abundance ra- 
tio (|Radzig fc Smirnovl . |1985) for 17 sub-components. We 
included 2 blends (Cri 4553.945 A and Zrii 4553.970 A) 
and 6 spectral lines of other elements observed in the far 
wings of this line. The upper panel of Fig. |2| shows that 
these blends produce only minor effects on the synthesized 
profile. 

We calculated the emergent intensities /aa along the 
Ban 4554 A line profile for the set of AA wavelength points 
and for every {xi,yi) vertical column of the 3D snapshot 
corresponding to the solar disc centre fi = cosQ = 1 (0 be- 
ing the heliocentric angle) . The profiles were normalized to 
the mean continuum intensity (Ic) averaged over the snap- 
shot. The continuum intensity was also used as a criterion 
to separate granular and intergranular regions. The {xi,yi)- 
grid points with continuum intensity greater than (Ic) were 
taken as granules (and the opposite for intergranules) . 

Since we are doing calculations for the disc centre, we as- 
sume complete frequency redistribution (CRD) for the Ban 
4554 A line. According to .RuttenI (|l978D . lRutten fc MilkevI 
(|l979t) the effects of partial frequency redistribution (PRD) 
increase towards the limb, where the frequency-dependent 
wing source structure becomes noticeable. At the disc cen- 
tre the differences between the PRD and CRD profiles of 
this line are small (1% of the continuous intensity in the 
core) so we can safely neglect the effects of PRD. 

We do not use any ad hoc parameters such as micro- 
or macroturbulcnce for the synthesis since the profiles are 
broadened in a natural way by the velocity field existing in 
the 3D model. The damping constant for the barium lines 
was determined as the sum 7 = 76 + 7rad of van der Waals 
collisional broadening 70 by neutral hydrogen and helium 
atoms and radiative broadening 7rad- The other collisional 
broadening processes ( Stark br o adeni ng, quadrupole broad- 
ening) are negligible (jRuttenl . 1X9781 ). We employ the 75 
based on a theory in which the van der Waals potential 
is replaced by a Smirnov-Roueff potent i al for the close in- 
teractions (jDeridder fc van Rensbergenl . Il976f ). In view of 
the uncertainty of collisional damping, we treat the 76-value 
as a free parameter by introducing the usual enhancement 
factor E. In this s tudy we use E = 1.3 a nd = 2.16 
derived earlier by lOlshevskv et al.l (|2008[ ) from the Ban 
4554 A 1.5D NLTE fine modelfing in the same 3D snap- 
shot. For these values of E and the authors obtained 
excellent agreement between the spatially averaged syn- 
thetic profile and the o bserved one taken from Liege atlas 
(|Delbouille et al.l . |1973( ) . We understand the limitations of 
the Smirnov-Roueff potential approximation for the close 
interactions. These limi tations might be overcome using the 
semi-classical theory of lBarklem fc O'Maral (j 19981) . The in- 
troduction of the broadening for P-D states given by these 
authors can be important for a proper modelling of the Ba n 
4554 A line, particularly, for any accurate calculation of the 
polarization Stokes amplitudes (Q/I) produced by scatter- 
ing processes in the solar atmosphere. The elastic collisions 
with neutral hydrogen atoms are known to be efficient in 
modifying the atomic level population of long-lived levels, 
like Ban 5^D. Recent investigation of the role of collisional 
depolariza tion of th e Ban 4554 A line in the low chromo- 
sphere by iDerouichI (j2008[ ) shows that this line is clearly 



through the effect of collisions on the 5^D level. However, 
the impact of depolarizing elastic collisions on the emer- 
gent intensity (/) profiles in a weekly anisotropic medium 
like the solar ph otosphere is expected to be negfig i ble (e .g., 
see the book by iLandi Degl'Innocenti fc Landolfil . |2004| ) . 

In Appendix [Al we summarize the NLTE effects on the 
formation of the Ban 4554 A line in the 3D snapshot and 
the compare NLTE and LTE calculations. 

Below, we use the concept of the Eddington-Barbier 
height of line formation, i.e. we evaluate the height H 
where the line optical depth at a given wavelength point 
A A is unity: t(AA) = 1. We then assume that the in- 
formation about vertical velocity Vz at wavelength point 
A A originates from height H. The reader should be 
aw are of the practical limita tions of such a concept (see 
San chez Almeida et al.l . Il996f) . However, for the purposes 
of statistical quantification of the errors in Doppler diag- 
nostic methods it is reasonable and convenient to use the 
concept of the Eddington-Barbier height of formation. We 
assumed that two-dimensional maps of the synthesized in- 
tensity profiles represent observations in the case of perfect 
seeing conditions and no instrumental effects. We then use 
such "perfect observations" to compute Dopplergrams and 
A-meter velocities and to compare them to the "true" snap- 
shot velocities at the corresponding heights H. 

2.3. Spatial smearing 

Since "perfect observations" are only of theoretical inter- 
est to an observer, we study here the case of several de- 
grading effects applied to the profiles. Observed images are 
degraded because of the Earth's atmospheric turbulence 
(seeing) and light diffraction b ythe telescope aperture (the 
finite spatial resolution of the telescope). There are other 
factors, such as stray light, but we expect that the observa- 
tions are corrected for them during the reduction process. 

Mathematically, the Fourier transform of the image reg- 
istered by detector is related to the Fourier transform of the 
original image via the modulation transfer function (MTF) 
as follows: 



/(/) = o(/) M/)), 



(1) 



where / is the spatial frequency in (rd~^), and () denotes 
the ensemble average. 0{f) is the Fourier transform of the 
observed object, /(/) is the Fourier transform of the im- 
age registered with a telescope camera, r(/ ) is the MTF 
of the atmosphere and telescope. Following iFriedl ()1966t ). 
for exposures longer than the characteristic lifetime of the 
atmospheric turbulence, the MTF is given by: 



{r{f))LE = ro(/)exp(-3.44(A//i?o)'/'), 



(2) 



where / = |/|, and ro(/) is the autocorrelation function of 
the telescope pupil: 



ro(/) = - 

TT 



COS 



D D \ \ D , 



1/2- 



(3) 



affected by isotropic collisions with neutral hydrogen atoms 



if Xf / D < 1. In such an interpretation only two param- 
eters define the resulting image quality: telescope diame- 
ter D and Fried's parameter Rq (Fried, 1966; Korff, 1973i). 
The latter parameter depends on seeing conditions and de- 
scribes characteristic size of atmospheric turbulence cells. 
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At each wavelength, the original 2D intensity maps were 
Fourier-transformed and multiplied by the MTF calculated 
for a known telescope diameter and Fried's parameter Rq. 
An inverse Fourier transform gives us the images registered 
by detector, i.e. the "observed" images, affected by the 
diffraction by the telescope aperture and the seeing effects. 



2.4. Determination of Fried's parameter 

In order to make a direct comparison between the synthetic 
and observed spectra we have to know Fried's parameter 
Rq. There is no direct way to me asure such a parameter in 
observations. iRicort et all (|l98 l) proposed determining Rq 



from the observed r.m.s. continuum contrast 6L 



of so- 



lar granulation. Their computations are based on using an 
analytical form for the power spectr um of the intens i ty dis - 
tribution at A = 5050 A obtained bv lRicort fc Aimd (|1979[ ) 
from speckle intcrfcrometric observations of the solar granu- 
lation. To establish calibrati on dependence betw een Rq and 
^-^r.m.s. at other wavelengths iRicort et al.l (|l98lh used sim- 
plified assumptions that the solar photosphere behaves as 
a blackbody, and that the differences in formation heights 
of continuum radiation as a function of wavelength can be 
ignored. 

Here we propose another way to set up calibration 
curves for determination of Rq from granular contrast in 
observations. Our calculations use a 3D radiative trans- 
fer solution for the solar continuum intensity obtained by 
iTrujillo Bueno k Shchukinal (|2009f ) in the same 3D snap- 
shot as was employed above for the Ban profile synthesis 
(see Sec. 12.11) . Note that such a 3D approach gives a possi- 
bility of avoiding the simphfication used bv IRicort et al.l 
(|l98lt) . The theoretical calibration curves for the wave- 
lengths between 4000 A and 8000 A are plotted in Figure [S] 

In order to verify the 3D approach we compare the em- 
pi rical calibrat i on cur ve at A = 5050 A taken from figure 1 
of lRicort et al.l (|1981h with our theoretical curve. The anal- 
ysis of the results presented in Fig. fallows us to conclude 
that at A = 5050 A both the empirical (open circles) and 
the theoretical (dash line) curves in fact coincide. Keeping 
in mind such a good agreement we use our theoretical 3D 
approach for determination of Rq in Section [3l 



2.5. Definition of artificial datasets 

Here we define the nomenclature for the artificial datasets 
to be used in the rest of the paper. 

DOT-LIKE DATA. ISuttcrlin et all (|2001l) used an 



19761 ) installed 3 



Irkutsk barium filter (jSkomorovskv et al. 
on the Dutch Open Telescope to obtain Dopplergrams in 
the Ban 4554 A line. We have calculated Dopplergrams 
from the synth e tic pr ofiles in order to model the results of 
IStitterlin et al.l ()2001f ). To do that, we smeared the orig- 
inal profiles using the DOT's diameter of £> = 48 cm. 
In addition, spectral smearing was performed by convolv- 
ing the profiles in wavelength with a Gaussian of FWHM 
(full width half maximum) of 80 mA, corresponding to the 
barium filter spectral profile. The Dopplergrams were con- 
structed as a second-order polynomial fit to the five points 
on the line profile situated at AA = —70, —35, 0, 3 5 and 
70 mA. The observations of ISiitterlin et all (|200lD were 
speckle-reconstructed, so we assumed that the influence of 





40 


: ' ' ' ' '/////// 


1 1 1' [ 1 1 / 1/ 


/' '//' /' '/' '/- 






- /////// / 




1 / / - 






; 


////// ' 


1 / / : 




30 




1/ / 


r / / / ~- 






' 1/ 


1 1 / 








: 8000 A li 


/////// / / 

/////'/ / A 


' / / Aooo A : 


B 




/ / / / 7 / // 




20 


~ //////// 

























X \ 




10 






\ 










D = 70 cm 











1 . . , , 1 . " 




5 


10 


15 20 








(51 (%) 

r.m.s. 





Fig. 3. Root-mean-square continuum contrast (5/r.m.s. as 
a function of Fried's parameter Rq for different wave- 
lengths and for telescope diameter Z? = 70 cm. (5/r.m.s. 
curves were calculated using the 3D radiative trans- 
fer formal solution for cont inuum intensity obtained by 
ITrujillo Bueno fc Shchukinal (2009) in the 3D HD model. 
Dashed and thick solid line: (5/r.m.s. at A = 5050 A and 
A = 4554 A, respectively. Thin solid lines: Wr.m.s. for 
wavelengths from A = 4000 A till A = 8000 A plot- 
ted with step 200 A. Filled circles correspond to Rq = 
10, 14, 19, 29, 38 cm. Op en circles represent the data for 
A = 5050 A obtained by IRicort et all (11981! ) using an em- 
pirical approach. 



the Earth's atmosphere was minimized, and that no smear- 
ing with the Fried parameter was applied. 

VTT-LIKE DATA. The second type of artificial data rep- 
resents observations done with the German Vacuum Tower 
Telescope (VTT) and Triple Et alon SOlar Spectrometer 
(TESOS; iTritschler erall . I2002D . The synthetic profiles 
were convolved using the telescope diameter D = 70 cm, 
and a spectral bandwidth of TESOS of FWHM=15 mA. 
The Fried parameters in this case ranged from Rq = 
10 cm (medium seeing) to Rq = 38 cm (excellent see- 
ing). According to Fig. [3] (thick solid line and filled cir- 
cles) the corresponding root-mean-square continuum con- 
trast (5/r.m.s. of solar granulation at the wavelength of the 
Ban 4554 A line varies from 7.5% to 16.5%. 



Comparison to observations 

In this paper we used observations of the Ban 4554 A 
line, obtained in September 2006 by E. Khomenko, M. 
Collados and R. Centeno at the 70-cm German Vacuum 
Tower Telescope (VTT) at the Observatorio del Teide in 
Tcnerife with the help of Triple Etalon SOlar Spectrometer 
fTESOS: l^itschler et aLl .l20021. The quiet Sun region close 
to the disc centre was observed, the presence of the mag- 
netic activity was controlled by the simultaneous obser- 
vation with the Tene rife Infrared Polarimeter II (TIP II; 
Collados et al.l . [20071 ) in Fei fines at 1.56 fim. We took a 
single TESOS wavelength scan in Ba 11 4554 A line made on 
September 2, which represents a series of 36 narrow-band 
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Fig. 4. Calculated and observed intensity disc-centre pro- 
files of the Ban 4554 A line. Solid line: the spatially aver- 
aged NLTE profile for the 3D model. Small open circles: the 
profile from the Liege Atlas (top); spatially averaged pro- 
file obtained from observations at the VTT with TESOS in 
2006 {bottom). 
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Fig. 5. Root-mean-square contrast at different wavelengths 
along the Ban 4554 line profile. Dashed line: Slr.m.sX^^) 
unsmeared synthetic data. Thick solid line: (5/r.m.s.(<5A) in 
the synthetic data smeared using Fried's parameter Rq = 
9 cm. Circles correspond to the observations at the VTT 
with TESOS in 2006. 



filter images (FWHM=15 mA) obtained with 10 mA in- 
terval along the 4554 A line profile. Each image was taken 
with a 500 /iS exposure, taking around 26 s to scan the 
complete line profile. The pixel resolution of the TESOS 
camera was 0.089 arcsec. Due to the instrument specifics, 



the TESOS field of view is circular. Taking into account 
that the image quality towards its edge is worse we have 
chosen a 256 x 256 pixel square cut from the central part 
of each image. By spatially averaging each of the 36 filter 
images we got the mean intensity profile of the Ball 4554 
A line norm alized to the mean con tinuum intensity of the 
Liege Atlas (|Delbouille et all . Il973f ). 

Figure m shows observed and computed spatially aver- 
aged profiles of the Ba ii 4554 A line. The top panel of Fig. |4] 
demonstrates that the profile computed using the 3D snap- 
shot is in a good agr eement with the observed one taken 
from the Liege atlas (jDelbouille et al.l . Il973l ) . The central 
part of the line is reproduced very well while the fit for 
the outer red wing of the computed profile is a little worse. 
The averaged line profile obtained in 2006 with the TESOS 
instrument at the VTT is shown in the bottom panel of 
Figured! Again, both computed and observed profiles agree 
well. The saw-tooth shape of far wings of the observed pro- 
file arises because of t uning effects of the TESOS etalons 
(|Tritschler et al.l . l2002f ). We would like to stress that such a 
good match of the observations is not based on free param- 
eter fitting but has been achieved self-consistently through 
the NLTE plane-parallel modeling combining a realistic 
atomic model and the 3D hydrodynamical snapshot. 

Figure O shows the run of the r.m.s. contrast with wave- 
length along the Ban 4554 A line for three sets of single- 
wavelength images: observed, smeared, and non-smeared 
synthetic data. The synthetic profiles were smeared employ- 
ing the method developed in Section [^31 We apphed the 
theoretical calibration curve at A = 4554 A shown in Fig. [3] 
(thick solid line) to define the Fried parameter from the 
VTT observations of the Ban 4554 A line. We found that 
during observations the r.m.s. continuum contrast (5/r.m.s. 
at the wavelength of this line was around 7%. Thus, we es- 
timate Rq to range range between 9 and 10 cm. Synthetic 
data in Fig. [5] are smeared using Rq = 9 cm. 

Seeing effects influence the r.m.s. contrast in the line 
wings, while in the line core its influence is not so signifi- 
cant, although the wing contrast is much higher than the 
contrast in the line core. Two local maxima near AA = ±50 
mA are not so high and sharp in the observed data as in 
the smeared computed proflles. While at other wavelengths 
the agreement is fairly good. Such blurring of the maxima 
may be due to the wider bandwidth o f TESOS fllter system 
than given bv lTritschler et all (|2002f ). 

Comparison of the synthetic and observed line proflles 
performe d in this section sugg ests that 3D atmospheric 
model of lAsplund et all ()2000af ) together with our NLTE 
calculations and treating of atmospheric and instrumental 
influence, reproduce well the average characteristics of the 
observed granulation pattern in Ba II 4554 A line. 

4. Dopplergram technique 

We discuss in this Section the range of validity of the 
5-poin t Dopplergram technique used by ISiitterlin et al.l 
(|200lD at the DOT to obtain LOS velocities. The compari- 
son between the Dopplergram velocities from the DOT-like 
data Vp and the "true" snapshot velocities Vz is given in 
Figure [6l The velocities Vp were corrected for the asymme- 
try caused by the hyperfine structure (the bottom panel of 
Figured]). 

When applying the Dopplergram technique, it is impor- 
tant to identify the height in the atmosphere from where 
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Fig. 7. Comparison of the snapshot vertical velocities Vz 
with the Dopplergram velocities Vp inferred from the 
"speckle-reconstructed" VTT-like data. 



the velocities originate. We assumed that the information 
about velocities at different points around the snapshot 
comes from nearly the same heights. Then we calculated 
the correlation coefficient between the maps of Vp and Vz, 
taking the latter at different heights, as shown in Figure[S^. 
This correlation tells us at what height the Dopplergram 
velocities of the DOT-like data measure the real solar ve- 
locities. As we can see, the correlation coefficient strongly 
depends on the atmospheric height, for both the smeared 
and non-smeared data. The maximum correlation reaches 
the value ~0.9 at heights around 300 km. Thus, we can 
identify a rather narrow layer where the Dopplergram ve- 
locities originate. 

Figure [SJa gives the scatter plot of Vp and Vz , the lat- 
ter taken at the height of maximum correlation (285 km). 
It demonstrates that in general the speckle-reconstructed 
Dopplergram velocities correspond to the true snapshot 
velocities rather well, with the standard deviation being 
less than 0.4 km/s. The surface maps of Vp and Vz also 
agree rather well (panels c and d). Thus, the Dopplergram 
velocities reproduce in many details the "true" granula- 
tion velocity structure existing at heights around 300 km. 
Nevertheless, it can be seen from the maps that at some 
locations in intergranular lanes the Dopplergram velocities 
overestimate the snapshot velocities. These are the loca- 
tions with strong downflows. The difference between Vp 
and Vz can be as high as 1 km/sec at these locations. Such 
an excess manifests itself in Dopplergram velocity maps as 
patches of enhanced brightness (Fig. [G}:, location marked 
by (1)), being absent in the snapshot velocity maps. 

Interestingly, with the narrower bandwidth and larger 
diameter of the telescope than in the case of the DOT obser- 
vations the overestimation of the Doppler velocities Vp in 
intergranular areas with strong downflows found from the 
"speckle-reconstructed observations" turns out to be appre- 
ciably larger. We display such a case for VTT-like data in 
Figure [T] As follows from the Figure, the deviation reaches 
approximately 3 km/s. 

Figure [5] explains why such bright patches can arise 
in speckle-reconstructed observations and why it happens 
specifically in intergranular areas. It displays examples of 



the parabolic fit to the five fixed wavelength points of the 
original profiles (top) and DOT-like smeared profiles (bot- 
tom) . The profiles are taken at the spatial locations marked 
(1) and (2) in Figure [Sb, c, d. The profile marked (1) orig- 
inates from the area with a high intergranular downflow 
(left) and the profile marked (2) originates from the area 
with a high granular upfiow (right). According to Fig. [5^, it 
seems impossible to measure correctly the velocity from the 
weak twisted intergranular profile provided that the strong 
redshift forces four of the five wavelength points used for 
fitting to be situated in the blue wing. As a consequence, 
the parabola turns out to be much more redshifted than the 
original profile. So the Dopplergram velocity derived from 
the wavelength position of the parabola bisector (thin verti- 
cal line) will be considerably greater compared to the veloc- 
ity inferred from the bisector of the original profile (small 
open circles). Smearing caused by the telescope and filter 
to some degree smoothes the irregular shape of the inter- 
granular profile and decreases the redshift of the parabola 
(Fig.[5)D). Nevertheless, the velocity remains overestimated. 
Thus, we conclude that the Dopplergram velocities calcu- 
lated from the speckle-reconstructed DOT-like data in in- 
tergranular lanes can be appreciably greater than the real 
ones. 

Contrary to the intergranular profile the deep and more 
symmetric granular profile allows us to estimate the ve- 
locity from the parabolic fit rather well (Fig. [8j:). The ef- 
fects of smearing are less pronounced in this case (Fig. [5]i) . 
So in granular areas with strong upflows one can expect a 
reasonable agreement between the "true" velocities Vz and 
velocities Vp recovered from the DOT-like data. 

Summarizing the conclusion of this Section, the five- 
point Dopplergram technique applied to the Ban 4554 A 
line profiles in speckle-reconstructed DOT-like data is a 
valuable tool for the diagnostic of the solar velocity field 
at heights around 300 km. Only in intergranular lanes with 
strong downflows can the velocity be overestimated produc- 
ing artificially bright points at the Dopplergram velocity 
maps like those obtained bv ISiitterlin et afl (j200lh . These 
authors interpreted the bright points as locations of the 
magnetic flux tubes. However, our calculations show that 
at least part of such bright points may be simply an artefact 
caused by the parabolic fit. 

5. A-meter technique 

The A-meter method of IStebbins fc Goodd (|1987D can be 
considered as a form of the "bisector shift" technique pro- 
posed earlier bv lKulander fc JefferiesI (fl966) to evaluate the 
atmospheric velocity field from the asymmetries in observed 
line profiles. The A-meter method deals with two line profile 
parameters such as the line bisector and the full spectral 
line width AAvy. 

The displacement of a midpoint of a section of the line 
profile of a certain width AXw is assumed to be a result 
of the Doppler shift of the line opacity coefficient caused 
by the non-thermal velocitie s in the layer where this sec- 
tion is formed. According to IStebbins fc Goodd (|1987 ^ the 
velocities and intensities from progressively deeper sections 
of the line profile correspond to progressively higher lay- 
ers in the atmosphere. Following the A-meter procedure we 
introduced a set of line profile widths AAvi/ ranging from 
306 to 78 mA (see Fig. [8]). Due to the irregular shape of 
the Ban 4554 A line core, a significant number of synthetic 
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Fig. 6. Comparison of the snapshot vertical velocities Vz with the Dopplergram velocities Vp inferred from the DOT-like 
data, (a): Correlation coefficients between Vp and 14 as a function of height for the case of no smearing (solid line) and 
telescope+filter smearing (dashed line), {b): Dopplergram velocities Vp vs. snapshot velocities Vz - The latter were taken 
at heights where the correlation between Vp and Vz reaches maximum. Blue crosses are granular points and red circles are 
intergranular points. The light blue curve is the average over bins with 100 surface points, (c): Map of the Dopplergram 
velocities Vp. Circles with numbers (1), (2) correspond to the spatial points where the emergent profiles have high red 
and blue shifts, respectively. Those profiles are shown in Figure [S] Red denotes downflows and blue, upflows. (d): Map 
of the snapshot velocities Vz ■ 



profiles never have a spectral line width lower than 78 niA. 
At the same time shallow far wings of many profiles do not 
allow us to detect and analyse velocities at widths above 
A\w ~ 306 mA. The selected A\w range corresponds to 
the intensities of spatially averaged profile varying between 
0.86 and 0.16. The corresponding range of the mean forma- 
tion heights lies between — 50 and ~ 500 km. 

We slightly modified the standard A-meter procedure 
keeping in mind that velocity field shifts the line opacity 
coefficients of various atmospheric layers towards one wing 
of the static line and away from the other. Such a shift 
expands the optical line depth r(AA) on one side of the 
line and compresses on the other side. As a consequence, the 
height of formation of the equal intensity points in the blue 
AAb and red AA^ wings belonging the same spectral width 
AXyy tends to be different. Our modified A-meter procedure 
is the following. We applied the above set of spectral widths 
both to the spatially averaged profile and to the individual 
profiles over the snapshot. The spatially averaged profile 
corrected for asymmetry due to granulation velocity field 
(see Fig. [21 bottom) was considered as a reference profile 
representing a stationary case. Since our study deals only 
with the relative changes of the Ba ii 4554 A line parameters 
no correction for the asymmetry caused by the hyperfine 



structure and isotopic shift was made. We derived the A- 
meter velocities Vxm in the blue wings by measuring the 
Doppler shift of the blue intensity point belonging a certain 
spectral width relatively the corresponding blue point of the 
spatially averaged profile (same applies to the red wings). 

The velocities Vxm obtained by this method were com- 
pared to the corresponding snapshot velocities Vz- The Vz 
values were specified separately for the blue AA;, and red 
AXr wavelengths belonging the same line width. These ve- 
locities were taken at heights where the optical depth at the 
corresponding wavelength points was equal to unity (see 
Sect. 2.2). The results of these calculation are displayed in 
Figs [S] to [H 

5.1. X-meter: original resolution 

Figure [5] shows the correlation coefficients between the 
maps of A-meter velocities Vxm and the snapshot velocities 
Vz taken at corresponding heights. It gives several cases for 
different seeing conditions (value of the parameter Rq) for 
the VTT-like data and also the case of no smearing. When 
the smearing is absent, the correlation coefficient for both 
blue and red wing velocities is rather high. It reaches a 
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Fig. 8. Ban 4554 A line profiles at two selected locations corresponding to a high downflow (left) and a high upflow 
(right), marked by (1) and (2) in Figur^ (a, c): profiles without smearing. Thin Hnes with dots: parabolic fit to the 
five wavelength positions AA = —70, —35, 0, 35 and 70 mA indicated by black circles. Open circles: bisector of the Ball 
profile; thin vertical line: bisector of the parabolic fit. The horizontal lines are full spectral line widths AXw — AA^ — AAf, 
specifying the wavelength distances between equal intensity points in the red and blue line wings, {b, d): same but for 
the DOT-like smeared profiles. Xi and are the locations of the selected grid-points in the 3D-snapshot. 
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Fig. 9. Correlation coefficients between A-meter velocities Vxm and snapshot velocities Vz as a function of the line 
profile width AXw- The snapshot velocities are taken at heights H where optical depth at the corresponding wavelength 
T(AAb(r)) equals unity (see Sect. 2.2). Each line corresponds to a different value of Rq. The uppermost solid line refers 
to the case of no smearing. Left panel: blue wing velocities; right panel: red wing velocities. The VTT-like data are used. 



maximum value of ^ 0.97 around AAvi/ — 200 mA, being 
slightly smaller in the inner wing of the line. 

The scatter plots and maps of the V\m and Vz veloc- 
ities in the inner wing of the line at AA^y — 78 mA are 
presented in Fig. [TU] for the case of no smearing. As follows 
from panel (a) of this figure, in the inner blue wing A-meter 
velocities agree essentially with the true snapshot velocities 
at that height in most of the grid points. The correlation 



coefficient between them is very high and is close to 0.9. 
The comparison of the velocity maps in panels (c) and (d) 
supports this conclusion. Nevertheless, there is a set of grid 
points belonging mainly to intergranular lanes where the fit 
is fairly bad. The map of the true velocities Vz in Fig. [TOH 
reveals the presence of a number of intergranular lanes (as 
defined by the continuum intensity) where the upflowing 
velocities are observed at heights of formation of the blue 
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Fig. 10. Comparison between the A-meter velocities V\m at AAu^ — 78 niA and snapshot velocities Vz, measured at 
corresponding heights H for the case of no smearing, (a): Scatter plot of Vxm and Vz in the blue wing. Black open circle 
marks the location Xi = 5.5 Mm, yi = 5.4 Mm. (6): Same for the red wing velocities, (c): Map of the blue wing V\m. (d): 
Map of the snapshot velocities Vz corresponding to the blue wing. The colour coding is the same as in Figure [H 



wing intensity at AA^^ = 78 mA. The A-meter velocities do 
not recover such upflowing points well, suggesting signifi- 
cantly smaller absolute value of the velocities or even their 
sign reversal. An example of such a location is marked by 
an open circle in Figure [TUl 

Figure [TUb demonstrates that in the red inner wing the 
match between Vz and V\m is appreciably worse. The A- 
meter velocities tend to be more redshifted. According to 
Fig. [nl the blue and the red wing V\m are essentially the 
same (except for a few grid points) , whereas the correspond- 
ing Vz are not. The reason for that lies in the different 
formation heights of the blue and red wing intensities at 
the same AXyy (we discuss this point in more details in 
Sect. 5.3). The velocities measured by the A-meter tech- 
nique correspond better to the heights where the blue-wing 
intensities are formed. 

Figure fT^ . b gives two more examples of the correlation 
between the A-meter velocities and the snapshot velocities 
in the outer wings of the line. The match is typically rather 
good and the correlation is high, except that the amplitudes 
of the A-meter velocities are systematically lower. The latter 
is easy to understand bearing in mind that the method gives 
the average information over a certain height range, thus 
leading to a decrease in the amplitude. 

In summary, under perfect conditions, the A-meter tech- 
nique allows us to obtain information about the LOS ve- 
locities with a rather good precision over the whole pho- 
tosphere, the agreement being a little worse for the inner 
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Fig. 11. Scatter plot of the red and blue wing A-meter ve- 
locities at AAn/ = 78 mA in the case of no smearing. Black 
crosses are granular points and light circles are intergran- 
ular points, the black open circle marks the point corre- 
sponding the location Xi — 5.5 Mm, yi — 5.4 Mm shown in 
Fig. [101 



wings of the line. The blue and the red wing V\m velocities 
are very close to each other. 
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Fig. 12. Scatter plot of the A- meter velocities V\m and the 
snapshot velocities Vz in the wings of the line at AXw = 
99 mA {left panels) and AXw = 196 mA (right panels) 
for the VTT-like data, (a and b): no smearing; (c and d): 
Fried's parameter Rq = 29; (e and f): Rq = 10. The gray 
scale coding is the same as in Fig. [11] 

5.2. X- meter: reduced resolution 

Apart from the case of the original resolution discussed 
above, Fig. [5] contains the calculation of the correlation 
coefficients between Vz and Vxm for the VTT-like data 
smeared to have a different spatial resolution by varying the 
Fried parameter from 38 to 10 cm. As expected, the corre- 
lation coefficient gets lower compared to the case of perfect 
seeing conditions. Nevertheless, even for the medium spa- 
tial resolution {Rq — 10 cm) the correlation for the inner 
wings {AXw = 100 mA) is rather high {^ 0.7). For the 
higher values of Rq the correlation coefficient increases up 
to 0.8-0.9. Opposite the case for no smearing, the correla- 
tion coefficient decreases from the inner to the outer wings. 
The velocity measured in the outer wings {AXw > 120 mA) 
is less precise. 

The effects produced by the spatial smearing on the 
amplitudes of the measured velocities AAvy are shown in 
Figure [T2l The following points can be underlined from this 
Figure: 



— In the original data the maximum absolute values of 
the velocities are higher above intergranular lanes than 
above granules. Spatial smearing reduces both the ab- 
solute values of the velocities and the asymmetry in the 
scatter plots. Nevertheless, the velocity asymmetry is 
still present. 

— Even under excellent seeing conditions {Rq = 29 cm) 
the outer wing A-meter velocities for AXw > 120 mA 
are less reliable that the inner wing velocities. 

— The inner wing velocities {AXw < 120 mA) still give a 
fairy good measure of the true velocities for the Fried 
parameter {Rq = 10 cm). 

— As the inner wings of the line are less sensitive to spatial 
smearing, the Ba ii 4554 A line can be useful to measure 
velocities mainly in the upper photosphere. 

5.3. Intensity formation heights along the line profile 

The A-meter technique only yields qualitative results as 
long as it is not accompanied by the knowledge of heights 
where the information on the velocity and intensity varia- 
tions comes from. In this section we give the results of the 
calculation of such heights. We calculated the Eddington- 
Barbier formation heights (see Sect. 2.2) for the intensity 
at each section of the Ba ii 4554 A line profile having a cer- 
tain spectral width AAvi/ . Below we use the notation Hi, for 
the blue wing intensity formation heights calculated at AAf, 
positions in the blue wing, and similarly for the red wing 
intensity formation heights H^, calculated at AA^ positions 
in the red wing. Both Ht, and Hr correspond to the same 
AXw- We repeated the calculation for the all grid points of 
the 3D snapshot. The results of such calculations are pre- 
sented in Figures [13] and [14] These figures give an answer 
to several important questions. 

Firstly, how realistic is to assume that Hb and Hr be- 
longing to the same AXw are constant over the 3D snap- 
shot? The results shown in Fig. [T3] suggest that this as- 
sumption is far from reality. At each fixed AXw, both Hi, 
and Hr vary in a rather wide range. In outermost sections 
of the profile {AXw = 306 mA) the heights vary between 
about —175 km and -1-75 km, while in the inner section 
{AXw = 78 mA) the range of the variations is larger and 
lies between -|-300 km and -1-800 km. Despite the large scat- 
ter, the mean formation heights of the each section of the 
profile (H) have a well pronounced dependence on AAw 
(thick solid curves in Fig. [T31) . This makes it possible to as- 
sign in, a certain way, a height dependence to the velocity 
measurements by A-meter technique. This can be done by 
ascribing the response of each AAvi^ section of the spectral 
line to the mean height (H). 

Secondly, is it correct to assume that the Hi, and Hr for- 
mation heights belonging to the same AAw are equal? In 
order to answer this question, we show in Fig. [14] two rep- 
resentative cases for Hi, and Hr calculated at two positions 
in the inner (AA^^ = 99 mA) and the outer (AAh^ = 196 
mA) wings of the line. The AXw = 196 mA corresponds to 
the section of the Ball 4554 A where the line opacity pro- 
file is shallow. So one can expect the difference in heights 
of formation Hi, and Hr caused by the Doppler shift of the 
opacity profile to be small. The width AXw = 99 mA be- 
longs to the inner section of the Ba ii 4554 A characterized 
by a steeper part of the line opacity profile. It is expected 
that the same Doppler shifts produce a more important dif- 
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Fig. 13. Formation heights of the intensity along the line profiles caleulated for a set of spectral widths in the blue (left) 
and red (right) wings for all grid points from the 3D snapshot. Thick solid lines and open circles: mean formation heights 
and their standard deviation, respectively. 



ference in the optical depths on opposite sides of the Ban 
4554 A line and, hence, a somewhat greater difference be- 
tween the Hb and Hr heights. What is the magnitude of 
this difference? 

As follows from Fig. [T4k. in the outer wings of the Ban 
4554 A line the blue Hb and red heights are indeed very 
close to each other, both varying over the 3D snapshot. 
The histograms displayed in Fig. [T4b demonstrate that the 
distributions of H}, and are similar and have a sharp 
cut off. It means that outer wings are formed in a narrow 
atmospheric layer with mean heights {Hb) and (i?r) of the 
blue and red wings being very close to each other. 

Figure fT4b shows a similar calculation for the inner sec- 
tions of the Ban 4554 A line. The Doppler shift of the line 
opacity profile leads to a complex behaviour of the Hi, and 
Hi, heights. Two groups of points can be distinguished in 
the scatter plot of Hr vs. H^ calculated over the snapshot. 
For the first group (symbols above diagonal), the H^ heights 
are concentrated in a rather narrow band around -1-200 km 
while the Hr heights extend for several hundreds of kilo- 
metres from nearly -1-100 km to -1-800 km. This group of 
points belong largely to the profiles coming from granular 
regions. For the second group of points (symbols below the 
diagonal) the behaviour is the opposite. This group origi- 
nates mostly in intergranular lanes. Only a relatively small 
number of the profiles have approximately the same H^ and 
Hr heights (symbols along the diagonal) . 

The histograms of the Hi, and Hr heights for the inner 
sections of the Ban 4554 A line (Fig.intl) are not symmet- 
ric. The histograms have a maximum at heights around 
-1-200 km and a sharp cut-off below this height. There is a 
long tail toward 600-800 km. This tail is more pronounced 
for the red wing heights Hr- Such an asymmetry suggests 
that the information in the inner wing comes from two 
distinct layers. In most of the grid points the intensity at 
AAvK = 99 mA is formed around -1-300 km, whereas in a 
smaller but still appreciable number of points the intensity 
comes from higher layers between -f 400 and -1-800 km. On 
the whole, the mean formation height of the inner red wing 
is larger than that of the inner blue wing. 



6. Conclusions 

In this paper we have analysed the range of valid- 
ity of the two Doppler diagnostic techniques using the 
Ban 4 554 A line, i. e . the 5-point Dopplergram tech- 
nique (ISiitterlin et al.l. 200ll ) and the A-meter technique 
(|Stebbins fc Goodd . ll987ir We have performed NLTE ra- 
diative transfer calculations of the Ban 4554 A intensity 
profiles in the 3D s napshot of hydrodynam ical simulations 
of solar convection (|Asplund et al.l . l2000al) , neglecting the 
effects of horizontal radiative transfer, but considering a 
realistic barium atomic model and taking into account the 
hyperfine structure and isotopic shift. The original resolu- 
tion profiles were smeared to reproduce the DOT-like and 
VTT-like data and study the effects of the limited resolu- 
tion into the reliability of the results produced by the two 
Doppler diagnostic techniques. 

Our results for the 5-point Dopplergram technique can 
be summarized as follows: 



— The NLTE simulations using the 3D hydrodynam- 
ical model support the opinion that the speckle- 
reconstructed Dopplergram velocities obtained from the 
DOT-like data give appropriate representation of the so- 
lar photospheric velocity field. 

— The information on the velocities obtained from such 
Dopplergrams comes from a thin atmospheric layer 
located at heights around 300 km. The speckle- 
reconstructed Dopplergram velocity maps reproduce in 
many details the "true" velocity structures existing in 
this layer. 

— The Dopplergram technique can overestimate the ve- 
locities in intergranular areas with strong downfiows. 
This excess appears in the speckle-reconstructed maps 
as localized points with enhanced brightness. The in- 
terpretation of such bright points in terms of magnetic 
fields has to be carried out with caution. At least some 
of them may be an artefact caused by the Dopplergram 
technique itself. 

Our conclusions regarding the A-meter technique are the 
following: 
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Fig. 14. (a, c): Formation heights of the blue wing intensities vs. formation heights of the red wing intensities at two 
spectral widths: AXyy = 196 mA and AXyy = 99 mA. Dark crosses - granular points; light circles - intergranular points. 
{b, d): Histograms of the formation heights of the blue wing intensities (solid line) and red wing intensities (dotted line) 
for the same AXw- 



— Under perfect seeing conditions the A-mcter technique 
allows us to obtain information about the LOS velocities 
throughout the photosphere with rather good precision. 
Only in the upper photosphere is the particular velocity 
structure with upflowing points in intergranular lanes 
not always well reproduced. 

— The velocities measured by the A-meter technique cor- 
respond better to the "true" snapshot velocities taken 
at heights of formation of the blue line wing, rather that 
the red wing. 

— Even for rather medium seeing conditions the inner 
wings of the Ba II 4554 A line give reliable information 
about the velocity field in the upper photosphere. The 
results from the outer wings arc less reliable. 

— The mean formation heights of each section of the Ba II 
4554 A line profile have a well-pronounced dependence 
on the spectral width of the section. This gives the pos- 
sibility to assign height dependence to the velocity mea- 
surements by the A-metcr technique. The non-thermal 
motions can be reliably measured with the A-meter tech- 
nique applied to the Ban 4554 A line throughout the 
photosphere up to the temperature minimum. 
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Appendix A: NLTE modeling 

Figure lA.ll (panels a 
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to f) shows the population depar- 
Bail A4554 A line source functions 
two spatial grid points of the 3D 
the typical granular and intergran- 
these models to illustrate the dif- 
results for granules and intergran- 
departure coefficients are defined as 



P = f^-NLTE / ?^LTE where tt-nlte and ulte are the NLTE and 
LTE atomic level populations, respectively. The Complex 
behaviour of the /3-coefBcients shown in Fig. lA.lb . b is 
a result of the interac tion of several NLTE mechanisms 
descri b ed in detail by Carlsson et "all (11993 ): iBruls et all 
(|1992[ ): IShchukina fc Trujillo Buenol (|200l[) Here we just 
point out that for the barium atom the most important of 
them are ultraviolet line pumping, ultraviolet overioniza- 
tion, resonance line scattering and photon losses. The reso- 
nance line scattering and photon losses manifest themselves 
as a divergence of the lower 6s^Si/2 and upper 6p^Pi/2, 
6p^P3/2 levels of the Ban resonance lines. This divergence 
results from the surface losses near the layer where the 
optical depth is equal to unity. The losses propagate by 
scattering to far below that layer. Interestingly, for the in- 
tegranule the divergence of the /3-coefficients arises in the 
innermost layers. This happens because the photon losses 
occur mainly through the line wings of the Ban 4554 A line. 
As follows from Fig. IA.2I the line wings in integranules are 
formed considerably deeper than in granules. Such a differ- 
ence in the formation heights is a result of the Doppler shift 
of the line opacity coefficient caused by the velocity fi eld. A s 
a consequence, in the intergranular model (see Fig. lA.U a, 
b) the divergence starts already in the lower photosphere 
while in the granular model it happens only in upper pho- 
tosphere at heights around 400 km. 

Another important conclusion that follows from 
Fig. IA.2I concerns the height of formation of the Ba n 4554 
A line. The lower departure coefficient /3; is close to unity. 
So the scaling of the line opacity with this coefficient can- 
not lead to an appreciable difference between the NLTE 
and LTE heights of formation of this line. 

The excess of Ba n ions at the levels with excitation po- 
tentials above 5 eV visible in the granule model is produced 
by the pumping via the ultraviolet Ball lines starting at 
6s ^S, 5(i^D, 6p^P levels. For the intergranule model the 
overpopulation arises only in the uppermost layers. Such 
behaviour of the /3-coefRcients corresponds to the tempera- 
ture stratification of the models. The overpopulation of the 
high excitation levels of Ban in granules occurs because 
here the excitation temperature of the ultraviolet pump- 
ing radiation field appreciably exceeds the electron tem- 
perature. In integranules such superthermal radiation, and 
hence the level overpopulation, is present only above the 
temperature minimum region. In addition, in the intergran- 
ular photospheric layers the photon losses in the ultraviolet 
lines are more pronounced than in granules. 

The Pu/ Pi ratio of the upper and lower level departure 
coefficients of the Ban 4554 A line sets the departure of 
its line source function Sl from the Planck function B. 
Figure rOK c. d) shows that this departure (reflecting the 
corresponding departure coefficient divergence in the upper 
panels of this Figure) is larger in the intergranular than in 
the granular model. 

Figure IA.3I demonstrates that such behaviour is typical 
also for the total source function 5tot at the wavelengths 
corresponding to the inner wings (AA < 76.5 mA). On av- 
erage, in intergranular regions it drops below the Planck 
function while in granules the effect is less pronounced. 
Moreover, in granular areas with strong upflows the total 
source function can exceed the Planck function. This ex- 
cess can be understood if one takes into account that the 
resonance source line function is described by the two-level 
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approximation, i.e. it approximately equals mean intensity 
J . In the regions with small photon losses (like granules) 
the J > B, hence, Sl and S'tot have to be greater than B 
as well. 

Figure ET](e, f) show the NLTE and LTE disc-centre 
line profiles for the individual granular and integranular 
models. The profiles displayed in Fig. lA.ll fg ) result from 
averaging of the emergent intensities corresponding only 
to the granular models. Averaged intergranular profiles are 
shown in Figure[XT](h). These two bottom panels quantify 
the statistical effect produced by the deviation from the 
LTE in two such types of the atmospheric models. The main 
conclusions that may be drawn from the results presented 
in Fig. lA.li IA.2i IA.3I are the following: 

— The source function deficit, as compared to the LTE 
assumption, is the main mechanism that controls the 
formation of the Ban A4554 A line. The line opacity 
deficit is small and, hence, unimportant. 

— The divergence between S'tot and B changes the shape 
of the individual profiles, particularly the intergranular 
ones. 

— On average, the deviations from the LTE lead to deep- 
ening (i.e. strengthening) of the spatially averaged Ban 
A4554 A line profiles. The NLTE effects are most pro- 
nounced around the line core and are generally more 
important in the intergranular regions than in the gran- 
ular ones. The mean difference between the NLTE and 
LTE line core residual intensities does not exceed 5% 
for granules and 10% for intergranules. 

— Towards the wings, the LTE becomes a valid description 
for the Ban A4554 A fine profile. 
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Fig. A.l. (a, b): Departure coefficients f3 of the Ban atomic 
levels vs. height in a representative granular (left) and 
intergranular (right) models. Thick solid and thick dash 
lines: the /3-coefficients for the ground 6s^Si/2 and upper 
6p^P3/2 levels of the Ban 4554 A line, respectively. Dash- 
dotted lines: the /3-coefficients for the upper level of the 
Ban 4934 A resonance line. Dotted lines: the /3-coefficients 
for the Ban levels with excitation potentials above 5 eV. 
(c, d): The line source function Sl of the Ban 4554 A line 
(dash) and Planck function B (solid) in units of the temper- 
ature vs. height for the same models, (e, f): The Ban 4554 
A line profiles for the representative granule and itergran- 
ule. Solid and dotted lines: NLTE and LTE, respectively, {g, 
h): Spatially averaged granular and intergranular profiles. 
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Fig. A.2. The NLTE (solid line) and LTE (open circles) 

heights of formation of the core of the Ban 4554 A line 
along the slice of the snapshot yi = 0.6 Mm. Dash-dotted 
line: continuum height of formation at 4554 A. Dotted line: 
height of formation of the line wing for the wavelength po- 
sition AA = —76.5 mA. The background image is the snap- 
shot vertical velocity V^. Negative (upflow) velocities Vz 
correspond to granules (dark), while positive (downflow) 
velocities to intergranules (light). 
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Fig. A. 3. The total source function S'tot(AA)/J3 measured 
in the Planck function units for three wavelength points 

AA = —76.5; —42.5; —8.5 mA situated in the blue wing of 
the Ban 4554 A line profile. Left: Maps of the Stot{^X)lB 
at the mean intensity formation heights at these wave- 
lengths. Horizontal line corresponds to the slice = 0.6 
Mm and the filled circles to the surface positions Xi = 1.8 
Mm (granule), Xi = 2.4 Mm (intergranule). Right: Scatter 
plots of Stoti^^) / B and velocities 14 for the same wave- 
lengths. The velocities are taken at heights of formation of 
continuum intensity. 



